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Abstract 

Background: Transposable elements (TEs) constitute an important source of genetic variability owing to their 
jumping and regulatory properties, and are considered to drive species evolution. Several factors that are able to 
induce TE transposition in genomes have been documented (for example environmental stress and inter- and 
intra-specific crosses) but in many instances the reasons for TE mobilisation have yet to be elucidated. Colonising 
populations constitute an ideal model for studying TE behaviour and distribution as they are exposed to different 
environmental and new demographic conditions. In this study, the distribution of two TEs, Osvoldo and Isis, was 
examined in two colonising populations of D. buzzatii from Australia. Comparing Osvaldo copy numbers between 
Australian and Old World (reported in previous studies) colonisations provides a valuable tool for elucidating the 
colonisation process and the effect of new conditions encountered by colonisers on TEs. 

Results: The chromosomal distributions of Osvaldo and Isis retrotransposons in two colonising populations of 
D. buzzatii from Australia revealed sites of high insertion frequency (>10%) and low frequency sites. Comparisons 
between Osvaldo insertion profiles in colonising populations from the Old World and Australia demonstrate a 
tendency towards a higher number of highly occupied sites with higher insertion frequency in the Old World than 
in Australian populations. Tests concerning selection against deleterious TE insertions indicate that Isis is more 
controlled by purifying selection than Osvaldo. The distribution of both elements on chromosomal arms follows a 
Poisson distribution and there are non-significant positive correlations between highly occupied sites and 
chromosomal inversions. 

Conclusions: The occupancy profile of Osvaldo and Isis retrotransposons is characterised by the existence of high 
and low insertion frequency sites in the populations. These results demonstrate that Australian D. buzzatii 
populations were subjected to a founder effect during the colonisation process. Moreover, there are more sites 
with high insertion frequency in the Old World colonisation than in the Australian colonisation, indicating a 
probable stronger bottleneck effect in Australia. The results suggest that selection does not seem to play a major 
role, compared to demography, in the distribution of transposable elements in the Australian populations. 



Background 

Transposable elements (TEs) are DNA sequences that can 
move along the genome and constitute an important frac- 
tion of the genomes of most organisms studied so far [1]. 
It is indisputable that genomes have evolved in close asso- 
ciation with TEs, which are characterised by high species 
specificity and high structural variability, even in closely 
related organisms. After several decades of experimental 
and theoretical work, the mechanisms controlling TE copy 
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numbers in populations and most aspects of their biology 
within species are still not well understood. Classical theo- 
retical models predict that autonomous elements can be 
present in the genomes of a population as long as their 
steady transposition rates can counterbalance their elimi- 
nation by natural selection [2,3] . This equilibrium hypoth- 
esis was considered the principal explanation for the 
maintenance of a constant TE copy number. However, not 
all populations are at equilibrium. For instance: first, epi- 
sodic transposition bursts promote TE genome invasion in 
natural populations exposed to different population and 
environmental regimes; second, in other population condi- 
tions such as weak selection [4] or migration rates lower 
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than transposition rates [5], it takes a substantial time per- 
iod to reach equilibrium; third, stress environmental con- 
ditions and other population phenomena, including 
demographic (as in our case) and mating processes, are 
among the prime causes for the distribution of TEs in gen- 
omes. Data concerning plants [6] and nematodes [7] 
demonstrate that demographic history and mating systems 
shape TE diversity. However, the exact way in which they 
contribute to TE genome distribution in Drosophila 
remains obscure, particularly because of the lack of com- 
parative studies on parallel populations submitted to simi- 
lar environment and/or population conditions. In 
particular, colonising populations subjected to genetic 
drift processes due to drastic decreases in population size 
(i.e. bottlenecks), or confronted by new ecological condi- 
tions, could constitute an excellent model for carrying out 
these studies and evaluate the effect of demography versus 
selection. 

Drosophila buzzatii, a cactophilic species that origi- 
nated in the Chaco region of NW Argentina, colonised 
the Old World and Australia 300 and 70 years ago, 
respectively. Whereas the former colonisation was prob- 
ably associated with South American colonial trade and 
is not well documented, highly documented information 
is available for the Australian colonisation. The intro- 
duction of D. buzzatii into Australia is associated with 
the program of biological control of the prickly pear 
(Opuntia sp.) infestation that devastated Australian agri- 
cultural crops in the 1920s [8]. D. buzzatii was most 
probably introduced into Australia inadvertently via rot- 
ten Opuntia cladodes containing larvae of Cactoblastis 
cactorurn, a natural moth parasite of Opuntia, largely 
being transported in shipments from Argentina, that 
was used for Opuntia pest control. Cactoblastis and the 
accompanying insects spread rapidly to all Opuntia- 
infested fields [8,9]. After the plague was controlled, 
Opuntia diminished drastically in population size, 
remaining established in scattered habitat patches [10]. 
Our previous studies concerning Old World D. buzzatii 
colonisation demonstrated a bimodal pattern of distribu- 
tion of the Osvaldo retrotransposon consisting of high 
and low insertion frequency sites. However, it is 
unknown whether this distribution pattern, valid for the 
Mediterranean colonisation that occurred three centu- 
ries ago, is valid for other colonisation processes. The 
aim of the present work is to study the genome distribu- 
tion of two retrotransposons, Osvaldo [11] and Isis [12], 
in some of these colonising Australian populations of 
D. buzzatii. Subsequently, these results will be compared 
to those of previous studies concerning the retrotranspo- 
son Osvaldo in colonising populations of the Iberian 
Peninsula and in the original Argentina populations in 
order to evaluate how far TE dynamics depend on popu- 
lation parameters inherent in colonisation processes. Our 



results show that Australian populations have a mixture 
of high and low insertion frequency sites for both TEs 
and the drift effect associated to colonisation is probably 
the main factor responsible for the TEs distribution fre- 
quency. High insertion frequency sites would represent 
insertions prior to the colonisation. Sites of low fre- 
quency would correspond to new insertions occurred 
after the colonisation. 

Results 

Genomic distribution of Osvaldo and Isis retrotransposons 

The distribution of Osvaldo and Isis was analysed in the 
whole genome of two D. buzzatii Australian natural 
populations. The shape of the chromosomal distribution 
of these TEs is depicted in Figure 1; Osvaldo and Isis dis- 
tributions are represented on chromomes 2 and 3, 
respectively, as an example. Both TE distributions are 
characterised by two kinds of sites: sites with low and 
sites with high (>10%) insertion frequencies. This pattern 
was not observed for Osvaldo on other chromosomes, 
where the frequency of the highest occupied sites never 
reached 10% (data not shown). In the case of Isis, the 
same pattern is observed on chromosomes 2 and 5, with 
occupancy frequencies up to 30% in some sites. Osvaldo 
has three (2B2a, 2D5a and 2G2f) and Isis five (2B2a, 
2E5a, 3B5f, 3E5b, 5Gla) high occupancy sites in the 
whole genome. The occupancy (number of times that a 
site is present in the population) profiles of the two ele- 
ments are similar in the two populations studied, but Isis 
has a slightly higher occupancy frequency than Osvaldo 
(Table 1). Table 2 summarizes the means and variances 
of copy number per haploid genome for the two retro- 
transposons. The mean copy number is low for both ele- 
ments in the two populations but higher for Isis than for 
Osvaldo. Overall mean copy number comparisons 
between chromosomes demonstrates that whereas the X 
chromosome presents the lowest mean copy number for 
both elements, chromosomes 2 and 3 harbor the highest 
mean numbers of Osvaldo and Isis copies, respectively. 
Deviations from Poisson distribution were tested using 
chi-square goodness of fit, pooling adjacent classes with 
low expected copy numbers. 

In order to test whether the TE insertions are distribu- 
ted in a repulsive (rarely inserted in nearby regions) or a 
contagious (frequently inserted in nearby regions) way in 
each chromosome separately, we computed the linkage 
disequilibrium for each pair of sites detected by analysing 
each possible 2x2 contingency table [13]. The observed 
distribution of disequilibria was compared with the 
expected distribution (assuming no linkage disequilibrium) 
using Fishers hypergeometric formula [14]. A battery of 
correlation tests was performed for each transposable ele- 
ment and chromosomal arm, and the values were pooled 
at intervals of 0.1. Figure 2 depicts the two tests that 
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Figure 1 Distribution of Osvaldo and /s/s on chromosomes 2 and 3 respectively. 
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produced significant values and represents observed and 
expected correlation distributions of Osvaldo and Isis on 
chromosome 3 of the Jandowae population. Tests were 
significant for the two elements owing to an overabun- 
dance of classes with positive correlation coefficients. In 
the case of Isis, an excess of the class -0.19;-0.10 is 
observed, but the test is within the limits of significance 
(P = 0.04) 



Comparative study of Osvaldo and Isis copy numbers on 
different chromosomes 

The factors controlling TEs in natural populations and 
the magnitude of their impact on host genomes are 
topics of much discussion. The genomic distribution of 
TEs has been interpreted as the result of negative selec- 
tion against their deleterious effects [3]. Negative selec- 
tion could act either against deleterious effects of 
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Table 1 Occupancy profiles of euchromatic sites in Australian colonising populations 

TE Populations Occupancy profiles 
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Population origin: Inglewood (IN) and Jandowae (JA). 

Occupancy profiles: number of times that each site is occupied in a population. For example, in the Inglewood population there are 37 sites occupied once, five 
sites twice, etc. 



insertions and/or of chromosomal rearrangements 
resulting from ectopic crossing over between TE copies 
[2], If selection plays a preponderant role in TE distribu- 
tion, fewer insertions should be present on the X chro- 
mosome than in autosomes because of the deleterious 
effects on the X chromosome in hemizygous males [3]. 
In order to investigate whether this negative selection 
affected the Osvaldo and Isis distribution in Australian 
populations, their copy numbers were compared 
between chromosome X and autosomes. Observed and 
expected proportions were compared using a G test [13] 
among autosomes (Gc), between X and autosomes (Gb) 
and among all chromosomes (Ga). It was assumed that 
the expected number of copies was proportional to the 
relative amount of DNA in each chromosome. The 
number of bands per chromosome in the cytological 
map of D. buzzatii [15] was considered directly propor- 
tional to the amount of DNA, as per previous studies 

Table 2 Test of the Poisson distribution of Osvaldo and 
Isis per chromosome and haploid genome 

Populations 
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TE 
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0.08 
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Isis 



x 
2 
3 
4 
5 

HG 



0,09 
0.75 
1.11 
0.19 
0.89 
3.02 



0.08 
0.75 
0.94 
0.18 
0.56 
2.66 



1.00 
2.56 
0.06 
4.13 
3.59 



0.07 
0.63 
1.23 
0.27 
0.62 
2.83 



0.06 
0.44 
0.98 
0.27 
0.46 
1.57 



2.92 
3.03 
0.26 

I. 31 

II. 90 



See Table 1 for population identification. TEs: transposable elements; Ch: 
chromosome; HG: haploid genome; m: mean copy number; V n : variance of 
copy number; df: degrees of freedom; Bonferroni's correction was applied to 
the data. *P < 0.05; **P < 0.01. The numbers of individuals analyzed are 
indicated in parenthesis. 



[16,17]. To estimate the expected TE proportions, the 
proportion of bands of each chromosome was multiplied 
by the total TE copy number detected in a population. 
The results of the G test, depicted in Table 3, demon- 
strate that differences in Osvaldo between chromosome 
X and autosomes are significant when sites of high 
insertion frequency are considered but not significant 
when these sites are omitted from the analysis. How- 
ever, Isis appears to be more strongly controlled by nat- 
ural selection than Osvaldo because the results of 
comparisons between chromosome X and autosomes 
were significant regardless of whether high insertion 
sites were considered. This element is characterised by 
lower and higher copy numbers on chromosomes X and 
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Figure 2 Observed and expected frequency distributions of 
correlation coefficients between all pairs of sites of 
chromosome 3 in Jandowae population. A) Osvaldo; B) Isis 
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Table 3 Comparisons of the proportion of Osvaldo and 
/sis among chromosomes 



Osvaldo 
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Chromosomes 


Exp. Prop. 
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JA 


IN 
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0.05 
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0.03 
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0.22 


3 


0.207 


0.15 
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34.67* 


*(9.37**) 
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125.28**(4.17) 


105.25* 
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3 


3.0 (2.36) 
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Exp. Prop.: Expected TE proportions; Ga: Comparison of the proportion of TEs 
among chromosomes; Gb: Comparison of the proportion of TEs between 
chromosome X and autosomes; Gc: Comparison of the proportion of TEs 
among autosomes; H: Heterogeneity; Df: degrees of freedom;*P < 0.05; **P < 
0.01; Bonferroni's correction was applied. Values of the test excluding high 
insertion frequency sites are indicated in parenthesis. 



3, respectively, than would be expected under a random 
genome distribution. A heterogeneity test [13] was car- 
ried out to determine whether the differences always 
occur in the same direction and the results demonstrate 
that the two populations generally do not differ for Isis 
and Osvaldo copy numbers over the chromosomes. 
Comparisons among autosomes, when the two popula- 
tions are pooled and sites of high insertion frequency 
are eliminated, gave significant results for the Isis ele- 
ment but not for Osvaldo. 

Inversion polymorphism in Australian populations of 
D. buzzatii 

The original South American populations of D. buzzatii 
present inversion polymorphisms on chromosomes 2 
and 4, with eight (2st, 2j, 2jz 3 , 2jq 7 , 2y 3 , 2jc 9 , 2r 9 , 2js 9 ) 
and two (4s and 4st) chromosome arrangements, respec- 
tively [18]. As a result of the colonising processes in 
Australia and the Old World, rare arrangements were 



lost and others have changed their frequency. The two 
Australian populations considered in this study present 
two arrangements in chromosome 2: st and j, which are 
also the most frequent in the original and the Old 
World populations. The frequencies of arrangements st 
and j are 0.40 and 0.60 in Inglewood and 0.36 and 0.64 
in Jandowae, respectively. In contrast to the Old World 
colonisation, arrangements 2jq 7 and 4s have never been 
detected in Australia [19]. Interestingly, the 2jz3 
arrangement, absent in our samples, was identified in 
other studies at a very low frequency [20]. 

Recombination has been considered the main factor 
responsible for the distribution of TEs on chromosomes 
[21,22]. According to the model proposed by Langley et al. 
[23], TEs can induce chromosomal rearrangements by 
ectopic crossing-over between homologous elements 
located in different chromosomal sites. If the model pre- 
dictions are true, we will find TE accumulation in regions 
of low recombination such as inside inversions and adja- 
cent to inversion break-points. In order to determine 
whether the arrangements are associated with high inser- 
tion sites for Osvaldo and Isis, we computed the product- 
moment correlation coefficients (r). The results depicted 
in table 4 demonstrate a significant negative correlation 
between several highly occupied sites of Osvaldo (2D5a 
and 2G2f) and the j arrangement in Jandowae. Interest- 
ingly, in the Inglewood population, correlation coefficients 
for Isis and Osvaldo have the same sign (except for the site 
2G2f of Osvaldo) for both elements, but neither is signifi- 
cant. Moreover, of the sites that have a significant correla- 
tion, only one is located inside the j inversion (2D5a); the 
other is located in its break point (2E5a). 

Comparisons of TE mean copy numbers between 
Australian and Old World colonisations of D. buzzatii 

The Old World colonisation by D. buzzatii has been 
studied intensively in our laboratory with the Osvaldo 
TE [16,17]. We combined these data with those from 



Table 4 Correlation coefficients between J chromosomal 
arrangement and high insertion frequency (HF) sites 

Populations 



Arrangement 



IN 



JA 



HF sites of Osvaldo 

2B2a 
2D5a 
2G2f 
HF sites of Isis 

2B2a 
2E5a 



Corr. 

0.02 
-0.15 
0.08 

-0.04 
0,09 



Corr. 

0.12 
-0.30** 
-0.29** 

-0,21 
0.02 



Corr: correlation coefficient. Only sites having a frequency >0.10 were 
considered. 



Garcia Guerreiro and Fontdevila BMC Evolutionary Biology 201 1, 11:1 1 1 
http://www.biomedcentral.eom/1471-2148/1 1/111 



Page 6 of 1 2 



the Australia colonisation to compare mean Osvaldo 
copy numbers between both colonising populations and 
the original ones. Table 5 summarizes the mean Osvaldo 
copy numbers in the natural populations from Argen- 
tina and the Iberian Peninsula studied previously, and 
those from Australia. We observed that the original 
Argentinean populations have the lowest mean copy 
number of Osvaldo, with a slight difference between the 
two sets of populations analysed in two previous works 
[16,17]. Colonising populations from the Iberian Penin- 
sula possessed the highest mean copy number. A Krus- 
kal and Wallis non-parametric test of analysis of 
variance was performed to investigate whether the mean 
copy number of Osvaldo differs among populations 
from a given area and between populations from differ- 
ent areas (Figure 3A). While no differences in copy 
number were observed among Australian populations 
(H = 0.53, P = 0.46) or among Iberian Peninsula popula- 
tions (H = 5.45, P = 0.36), there were significant differ- 
ences among populations from Argentina (H = 16.17*, 
P = 0.02). When the three population groups were com- 
pared, the differences were significant (H = 219.84**, 
P = 0.00) and the Australian populations had mean 
copy number values intermediate between the Argenti- 
nean and Old World populations. The same analysis 
was carried out excluding positions of high insertion fre- 
quency in the colonising populations (see Figure 3B); 
the significance disappeared in all cases except when 
Argentinean and Old World populations were compared 
(H = 7.27**, P = 0.007). We conclude that the high 
mean copy number in colonising Australian populations 
can only be explained by the presence of some sites 
with high insertion frequency. Colonising populations 
could have suffered a genomic redistribution of the 
Osvaldo retrotransposon, probably due to the founder 
effect. As a consequence, certain copies arriving with 
the founders have a high occupancy in the current 
populations. These results agree with the hypothesis 
that the mean copy number of Osvaldo increases during 
the colonisation process. However, in the case of Old 
World colonising populations, differences still persist 



after elimination of high insertion frequency sites. The 
most likely explanation for this is the large mean copy 
number heterogeneity among Argentinean populations, 
because when the two populations with the lowest copy 
numbers are eliminated from the comparison, the differ- 
ences become insignificant. In the case of Isis, a recently 
characterised element, no information is available from 
previous population studies, but the mean copy number 
of laboratory lines from Spanish populations is approxi- 
mately three copies per genome [12], a value similar to 
that found in Australian populations. 

The sites with the highest insertion frequency of 
Osvaldo in different colonising populations from the 
Iberian Peninsula and Australia are presented in Table 
6. There was a general tendency towards lower numbers 
of highly occupied Osvaldo sites in Australia than in the 
Iberian Peninsula. 

Discussion 

TE copy numbers and distribution: Australia vs. Old World 
colonisations 

Comparing this study with previous works that analysed 
original (from Argentina) and colonising (from the Iber- 
ian Peninsula and Australia) populations, some differ- 
ences are apparent. While the original populations 
always demonstrate low occupancy per chromosomal 
site, the colonising populations (Australian and Spanish) 
have some sites with high occupancy. When colonising 
populations are compared, Old World populations have 
more highly occupied Osvaldo sites (eight) and a higher 
occupancy rate of some sites (up to 63% occupancy) 
than the Australian populations (three sites with less 
than 33% occupancy). The mean copy number is lower 
in Australian than in Old World populations, placing 
Australia in an intermediate position between Argentina 
and the Old World. Data from Isis are limited to the 
Australian populations, but interestingly, the mean copy 
number is similar in Australian and Old World labora- 
tory stocks. 

The overall picture that emerges from these data sug- 
gests a global increase in the mean copy number per 



Table 5 Mean copy number of Osvaldo retrotransposon in different natural D. buzzatii populations: original and 
colonising 



ELEMENT 




POPULATIONS 


Osvaldo 


Original 


Colonising 



Argentina a Argentina b I. Peninsula a I. Peninsula b Australia c 

(6) (2) (5) (1) (2) 



Mean 1.22 0.85 2.85 2.99 1.75 

The numbers of populations analyzed are in parenthesis; subscripts represent a detailed explanation about populations and bibliographic references; a : 
Argentinean populations (Amaicha, Arroyo Escobar, Catamarca, Cebila, Palo Labrado and Vipos) and Iberian populations (Carboneras, Sanlucar, La linea, Mazarron 
and Albufeira) as described in Labrador et al. 1998; b : Argentinean populations of Berna and Puerto Tirol and Iberian population of Carboneras, as described in 
Garcia Guerreiro and Fontdevila; c : Australian populations (Inglewood and Jandowae) presented in this work. 
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Figure 3 Kruskal-Wallis non parametric test of analysis of variance for the mean copy number of Osvaldo. A) All sites included, B) 
Excluding high insertion frequency sites. Populations identification: JA (Jandowae), IN(lnglewood), Argentinean and Spanish populations are 
explained in references [16,17]. 



site of Osvaldo during the colonisation process. This 
may be explained either by the drift effect following the 
founder event acting on many low occupied sites (turn- 
ing them into high occupied sites) or by increases of 
site-directed transposition rates in colonisation, or both. 
According to previous work [16,17,24,25], sites 



demonstrating high insertion frequencies are most likely 
due to a founder effect that occurred during the coloni- 
sation process, and low insertion sites are probably the 
result of new transposition events. 

The Australian and Spanish populations have similar 
insertion-site distribution patterns but some of the 
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Table 6 The highest frequencies of Osvaldo insertion 
sites in two groups of colonising populations 


High frequency sites 


1. Penmsula a (6) 


Australia b (2) 


2B2a 


0.48 


0.31 


2D5a 




0.18 


2D1h 


0.22 




2F4a 


0.42 




2F1fg 


0.18 




2G2f 




0.32 


2G2h 


0.26 




3F2b 


0.16 




4C2g 


0.10 




5A4b 


0.63 





The numbers of populations considered are in parenthesis; subscripts 
represent bibliographic references: a Labrador et al. 1998; Garcia Guerreiro and 
Fontdevila, c this work. See Table 5 for population identification. Only sites 
having a frequency >0.10 are presented. 



aforementioned differences deserve more detailed analy- 
sis. First, high insertion frequency sites are generally 
located in different chromosomal positions in the two 
colonisations, indicating two independent colonisation 
events from Argentina. Second, the number of highly 
occupied sites is greater in the Spanish than in the Aus- 
tralian colonising populations. These differences could 
be related to specific differences in the two colonising 
processes. While Old World colonisation appears to be 
caused by a low number of founders [17], the number 
of D. buzzatii individuals introduced during the Austra- 
lian colonisation was certainly large initially [8,9] but fell 
dramatically following the Opuntia biological control 
program. Therefore, we can imagine a scenario where 
D. buzzatii spread rapidly through the Opuntia area 
until 200,000 hectares of prickley-pear had been largely 
destroyed by biological control [26]. Since D. buzzatii is 
restricted to the Opuntia distribution area [27], a paral- 
lel drastic decrease of the D. buzzatii population size 
should have occurred. This rapid reduction in the num- 
ber of coloniser survivors could explain, at least in part, 
why the number of highly occupied sites in Osvaldo is 
lower in the Australian colonisation than in the Old 
World one. 

The results showed here are in accordance with pre- 
vious comparisons using allozyme variability [19,28], 
microsatellites [10,29] and mitochondrial DNA [30], 
where a decrease in variation was observed in Australian 
populations, probably due to a bottleneck followed by a 
population expansion after the initial founder event. 
Though the founder effect hypothesis is widely accepted 
by all authors working on Australian D. buzzatii popula- 
tions, there are discrepancies in relation to the bottleneck 
magnitude depending on the markers used. Whereas 
nuclear markers such as microsatellites [10] and allo- 
zymes [31] suggest a moderate bottleneck, mitochondrial 



DNA [30] data point to a very low number of founders. 
The mtDNA was considered a particularly valuable mar- 
ker for studies of recent population history owing to its 
maternal inheritance, lack of recombination and rapid 
evolution (relative to nuclear genes) [32,33]. The fact that 
data concerning TE distribution tend to be more in 
agreement with mitochondrial data, and that records of 
the cactus control program implied a strict quarantine of 
associated insects, suggest that the most likely explana- 
tion is that only a few individuals founded the Australian 
populations. Therefore, it is expected that the magnitude 
of the founder effect was stronger in Australian popula- 
tions than in Old World populations 

On the other hand, during colonisation selection could 
also play un important role for the elimination of deleter- 
ious TE insertions. For example, homozygosity can 
increase during the colonisation process, if a few coloni- 
sers are involved in the process, increasing the strength 
of selection against deleterious insertions as observed in 
selfing plant species [34,35] Additionally, selection could 
contribute to the adaptation to new environments as 
seen during the expansion of D. melanogaster from 
Africa [36,37]. Recent results demonstrate a role for 
some TEs in the adaptation of D. melanogaster colonising 
populations to the out of Africa temperate climate [38]. It 
is not impossible that some insertions are maintained 
because of their positive effects in terms of adaptation to 
the new environments encountered by colonisers. The 
bottleneck suffered by Australian populations could have 
produced gene frequency shifts that have been proposed 
as responsible for driving populations to new adaptive 
peaks by selection. But, in the present study, the presence 
of high insertion sites is more likely to be due to the drift 
effect associated to the colonisation process because the 
Australian D. buzzatii colonisation time of approximately 
840 generations (70 years multiplied by 12 generations 
per year) is probably insufficient to detect selection 
effects unequivocally. This hypothesis is supported by 
results in previous studies where the Osvaldo DNA 
altered sequence and its flanking genomic sequences in 
high frequency sites were identical in all Old World colo- 
nising populations [24], a result highly improbable under 
the selective hypothesis Thus, when considering all 
results together we can conclude that selection does not 
seem to play a major role, compared to demography, in 
the distribution of transposable elements in the Austra- 
lian colonising populations. 

TE dynamics in the colonisation processes 

Colonisation has been suggested as one of the processes 
responsible for awakening transposable elements in 
D. simulans [39]. During colonisation, individuals are 
exposed to stress attributable to new environmental condi- 
tions and population regimes. Environmental stresses can 
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induce changes in chromatin structure [1] and promote 
TE mobilisation as observed in plants [40] and other 
organisms [41]. D. buzzatii, a cactophilic species, has spe- 
cific resource requirements that imply intense natural 
selection for survival in extreme environments. Therefore, 
these species are less able to avoid environmental stress by 
habitat selection [42]. 

Considering a very simple model of colonisation, trans- 
position rates in colonisers can be estimated. We assume 
that sites originating from Argentina during colonisation 
are found at high frequency in colonising populations and 
that new transpositions correspond to low frequency sites. 
The estimations begin with the historical data, which sug- 
gest that Australian colonisation of D. buzzatii occurred 
approximately 70 years ago [31,43] and assume a mean of 
12 generations per year. Taking the pooled data from the 
two Australian populations (155 genomes), we estimate 
840 generations of colonisation (70 years multiplied by 12 
generations), 93 new sites and a mean copy number for 
Osvaldo insertions of 1.75. We estimate a minimum value 
of transposition rate per site per generation in colonisers, 
for Osvaldo, of 4.08 x 10" 4 (93/155 x 1/840 x 1.75) and for 
Isis of 3.62 x 10" 4 . These values have the same order of 
magnitude as those calculated for Osvaldo (1.28 x 10" 4 ) in 
previous studies with natural populations from the Old 
World [17]. No differences in transposition rates were 
observed between the two colonisations, or between the 
two transposable elements considered in this study. 

We cannot completely disregard local transposition 
events in some coloniser populations. It is well known 
that transpositional bursts have occasionally been found 
in Drosophila melanogaster laboratory stocks [44], other 
Drosophila species [45] and in plants subjected to differ- 
ent environmental stresses [46]. However, transposition 
bursts are difficult to see in nature in real time because 
TE copies involved in these events are later silenced 
owing to their deleterious effect to the host. The major- 
ity of these mechanisms include host factors as methyla- 
tion [47], deletions [48], silencing by RNA interference 
[49] or random mutations, excisions and purifying selec- 
tion [50]. 

In the case of D. buzzatii colonisation processes, we can 
accept that the populations analysed have not suffered 
large transposition rate changes. Another possibility is to 
imagine that even if colonisation induced local transposi- 
tion it could be unnoticed owing to the time elapsed, espe- 
cially if the TE eliminations are quick, as observed in the 
Helena element in genomes of D. melanogaster and 
D. simulans [48]. Old World colonising populations of 
D. buzzatii contain TEs where active and inactive copies 
coexist [24] and the low occupancy sites observed could 
be the result of a unique transposition or very few master 
copies present in the genome. However, the possibility 
that some low occupied sites harbored copies that have 



diminished their frequency in populations because of drift 
or selection effects against insertions with a deleterious 
effect cannot be disregarded. A colonising population 
from the Old World (Carboneras) demonstrated a 
decrease in the insertion frequency of many highly occu- 
pied sites (2B2a, 2F4a). In other studies concerning D. sub- 
obscura colonising species, the existence of common low 
frequency sites of bilbo and gypsy elements was attributed 
to a decrease in insertion frequency in some sites in colo- 
niser populations [25]. However, we hypothesize that high 
occupancy sites correspond to sites that have increased 
their frequency by drift associated to the colonisation pro- 
cess. Arguments in favor of this hypothesis are that sites 
of both elements, found at a low frequency in the original 
populations, are now at a high frequency in populations 
resulting from the two colonisations. However, a common 
highly occupied site (2B2a) of Osvaldo and Isis elements 
was observed in the Australian populations. This site 
could be a transposition hot spot induced by the colonisa- 
tion process as it is present at a very low frequency in 
some original populations. Moreover, the two elements 
are inserted in the same chromosomal band, confirming 
previous results where Isis retrotransposon appeared to 
have preferential insertions inside Osvaldo sequences [12]. 

Comparisons of TE copy numbers between chromosome 
X and autosomes, after elimination of high insertion sites, 
demonstrated that Isis appears to be controlled by purify- 
ing selection in Australian populations, whereas the selec- 
tion effect is not so evident for Osvaldo. However, the 
comparison is highly significant when the two Australian 
populations are pooled, which may suggest a stronger 
selection intensity against Osvaldo insertions in these 
populations. These differences may be due to a stronger 
selection effect in the Australian colonisation than the Old 
World one. Alternatively, the differences observed between 
the two elements could be due to the existence of differen- 
tial transposition and regulation mechanisms. Another 
explanation may be that selection effects on the Osvaldo 
element could go unnoticed if transposition events 
occurred only a few generations ago. A third possibility is 
that selection acts more strongly on the Isis element (for 
example, if the element is inserted inside genes or regula- 
tory regions) than on Osvaldo and, provided the Australian 
colonisation is recent, its effects are only detected in the 
element under the most selection pressure. 

Comparisons of the proportion of elements between 
autosomes demonstrated large differences due to an over- 
representation of Osvaldo on chromosome 2 (58-63% 
from the total) and of Isis on chromosome 3 (37-44% from 
the total). D. buzzatii Australian populations have one 
inversion on chromosome 2. Inversions [23,51], inversion 
break-points and nearby regions [52-54] are considered by 
some authors as sites of accumulation of TEs. These 
regions are characterised by a reduction in recombination 
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rates that diminishes the probability of deleterious ectopic 
exchanges. Correlations between highly occupied sites of 
Osvaldo and the J arrangement produce negative values; 
only two sites were significant: one located outside and 
another inside the J inversion. The most likely explanation 
is that most correlations observed are a consequence of 
the founder effect. This would explain why significant 
negative correlations are observed even in sites located 
inside the inversion. A similar interpretation was evident 
in a study concerning colonising D. subobscura popula- 
tions, where the founder hypothesis was favored by the 
fact that all correlations between sites and arrangements 
were significant only in the colonising populations [25]. 
Moreover, there were positive associations between chro- 
mosomal arrangements and highly occupied sites located 
outside of inversions. This result is not unexpected if we 
take into account that the effective population size is 
reduced during the first stages of colonisation. If the TE 
abundance is regulated by the effect of ectopic recombina- 
tion between elements [2], the existence of insertions in a 
homozygous state could reduce ectopic exchange events 
[35,55] outside inversions, leading to a relaxation of selec- 
tion pressure. 

Conclusions 

Knowledge concerning TE distribution in natural 
populations, subjected to different environmental and 
demographic regimes, is an important step towards 
understanding their role in species evolution. Here we 
demonstrate that the two retrotransposons studied in 
Australian populations presented bimodal distribution 
with high and low insertion frequency sites. These 
results indicate that Australian D. buzzatii populations 
were subjected to a founder effect and a strong bottle- 
neck during the colonisation process. This work 
emphasizes the importance of population factors in 
remodelling the distribution of TEs in the genomes of 
natural populations. 

Methods 

Drosophila strains and mating system 

Two natural populations from Southeast Australia: Ingle- 
wood (31.55° S/151.52°E) and Jandowae (26.46° S/151.50 0 
E), were caught in March 2006. Of the individuals col- 
lected, approximately one hundred males and 15 females 
of D. buzzatii were checked. Immediately after their arrival 
at the laboratory each wild male was crossed with three 
virgin females from the control line 63 42/7 F81 [24]. Iso- 
female lines were established from each wild female and 
subsequently one Fl male crossed in the same way as 
described for wild males. The control line was checked by 
in situ hybridisation before the experiment and was char- 
acterised as devoid of euchromatic copies of Osvaldo 



bearing two copies of Is is element (2Fld and 4F2b), with 
constant insertion profile over different generations. 

Only female larvae were analysed in crosses with the 
control line in order to have the whole wild haploid gen- 
ome, including X chromosome. The Is is insertion profile 
of each individual was obtained by subtracting the inser- 
tion profile of the control line from that of the Fl larvae. 
The minor disadvantage of using this method is that, in 
the case of Isis, the putative common sites between wild 
individuals and the control line were omitted from our 
studies. However, this is a powerful and reliable method 
used in previous studies concerning transposable element 
distribution inD. buzzatii andD. subobscura [16,17,25] 

DNA probes 

A 6.4 kb fragment of the Osvaldo retrotransposon, 
including the main part of the complete element except 
the two LTRs, was used and cloned in the pBSK vector. 
The Isis probe consisted of a 3.4 kb PCR fragment 
including pol-env [12] inserted in the pGEM-T vector 
(Promega). Both probes were labeled with digoxigenin 
lld-UTP (Roche) using a random primer reaction. 

In situ hybridisation 

Polytene chromosome squashes were prepared from D. 
buzzatii salivary glands of third-instar larvae as 
described by Labrador et al. [56]. Prehybridisation, 
hybridisation and post-hybridisation washes were per- 
formed following a protocol by Roche [57]. 

In situ hybridisation is a reliable technique and the most 
suitable method for detecting and locating insertions in 
chromosomal arms. However, the power of resolution of 
this technique is limited in cases of discrimination of very 
close sites on the chromosome or elements whose sequence 
has a divergence below 10%. In the present study, where 
only a few sites were detected, this is not a drawback. 

Inversions 

Inversions were identified after microscopic observation 
of the same slides used for in situ hybridisation. Inver- 
sion frequencies were calculated by dividing the number 
of inversions of each type by the total number of indivi- 
duals analysed. 

Statistical analyses 

Most statistical analyses and graphs were performed 
using the statistical software SPSS version 15.0 and 
Excel. In cases of multiple testing, Bonferroni's correc- 
tion was applied to the data [58]. 
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